DMPK Equation for the Edge Transport of Quantum Spin Hall Insulator 
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Using the random matrix theory, we investigate the ensemble statistics of edge transport of 
a quantum spin Hall insulator with multiple edge states in the presence of quenched disorder. 
Dorokhov-Mello-Pereyra-Kumar equation applicable for such a system is established. It is found 
that a two-dimensional quantum spin Hall insulator is effectively a new type of one-dimensional (ID) 
quantum conductor with the different ensemble statistics from that of the ordinary ID quantum 
conductor or the insulator with an even number of Kramers edge pairs. The ensemble statistics 
provides a physical manifestation of the ^-classification for the time-reversal invariant insulators. 

PACS numbers: 73.63.Nm, 73.61.Ng, 72.10.Bg 
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One of the recent advances in condensed matter 
physics is the discovery of the quantum spin Hall in- 
sulator (QSHI)i2. QSHI is a new type of topological 
insulator, which is gaped in the bulk, but has gapless 
edge modes that give rise to a quantized conductance. 
The key theoretical observation is the ^-classification 
for the time-reversal (TR) invariant insulating systems^: 
a two-dimensional (2D) insulator with an odd number of 
Kramers pairs of edge states and that with an even num- 
ber are topologically distinct, and the QSHI has an odd 
number of Kramers pairs at its edge. Such a classifica- 
tion has been established by the analyses on the topolog- 
ical structure of the Bloch bands 3,4 , and its robustness 
against the imperfections, such as the electron-electron 
interaction^ and disorders^^, has also been discussed. 
Experimentally, a quantized conductance is observed in 
HgTe quantum wells, and is taken as the signature of 
the QSHI phased, albeit not conclusively. Other exper- 
imental techniques, such as ARPES, are also employed 
for searching the new QSHIs^. At present, it is highly 
desirable to have more associations between the abstract 
^-classification and the physically measurable proper- 
ties. 

In this Rapid Communication, we investigate the en- 
semble statistics of the edge transport of QSHI in the 
presence of quenched disorder. In essence, a two- 
dimensional (2D) QSHI is effectively a one-dimensional 
(ID) quantum conductor with an odd number of Kramers 
pairs of conducting channels. Such a ID quantum con- 
ductor is actually a new species that can only be real- 
ized at the edge of a 2D QSHI^, different from the or- 
dinary ID conductors which always have an even num- 
ber of Kramers pairs of conducting channels. We estab- 
lish the Dorokhov-Mello-Pereyra-Kumar (DMPK) equa- 
tion 1 ^ applicable for such a system, upon which the en- 
semble statistics of the edge transport of the QSHI is 
investigated. The distinct ensemble statistics of the edge 
transport of the QSHI presents a physical manifestation 
of the .^-classification, and could be a useful probe for 
identifying the new TR invariant topological insulators. 

We consider a configuration shown in Fig. [Uja). Be- 
cause the insulating bulk prevents the direct communi- 
cation between the two edges, the system can be con- 




Figure 1: (color online) (a) The geometry of the system: an 
insulator with multiple edge conducting channels (represented 
by the color shaded bands at the two edges) is connected 
to the left and right measurement leads, (b) The transport 
along each of edges can be characterized by a S-matrix. af 
(b+) and a~ (b~) denote the right-going and left-going wave 
amplitudes, respectively, a* (or b i ) with the same index i 
are related by the TR and form a Kramers pair. TV denotes 
the total number of Kramers pairs at each edge, and can be 
odd (QSHI) or even (ordinary insulator). 



sidered effectively as two independent ID quantum con- 
ductors arranged in parallel. Each ID conductor has N 
Kramers pairs of conducting channels. We assume that 
the spin-orbit coupling is present, so the spins are not 
conserved in general. We do not assume the origin of 
the edge modes: they can be a result of the topological 
structure of the bulk bands, or from the extrinsic origins 
such as the surface dangling bonds. 

In general, the transmission along the ID conductor 
can be characterized by a 2N x 2N 5-matrix, which re- 
lates the incoming (ipm) and outgoing ('(/'out) wave am- 
plitudes: 
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,b^) T (see Fig.QJb)). In our la- 



beling of the channel numbers, TR symmetry imposes 
the constraint on the /S-matrix-i 2 -: 

S T = -S, (2) 

Moreover, the current conserving implies /S-matrix must 
be unitary: S^S = I. 

Under these constraints, the polar decomposition of 
the S-matrix read o 11 ' 12 : 



(3) 










£ 


r " 




' u 


" 





V T _ 




-T 


-s 







V 



2 



where U and V are TV x TV 
is a block diagonal matrix £ 
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diag[V7T, vTT, 
N = 2n (odd TV = 2n + 1), Tj denotes the ith transmis- 
sion eigenvalue. One immediately sees that for the odd 
TV, there is always one conducting channel that has the 
perfect transmission, without being adversely affected by 
the disorder. This is the reason behind the robust edge 
transport of QSHI. 

The different ensemble statistics for the ID conduc- 
tors with the odd and the even TV can already be ob- 
served if we compute the volume element of the con- 
figuration space expanded by the independent parame- 
ters of the S'-matrix: d/j,(S) — Jodfx(U)dfx{V) Yii dT,. 
To get the invariant measure Jo, we c alculate ds 2 = 
Tr[dS' t dS'] = g^dxidxj, and Jo = y / det[<^], where 

13 



= {U l] ,V l3 ,T i y^. We obtain 
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To determine the ensemble statistics of the ID con- 
ductor, we derive the Dorokhov-Mello-Pereyra-Kumar 
(DMPK) equation for evolution of the joint probabil- 
ity distribution function of the transmission eigenvalues 
about the length L: P(Ti,T2 ■ ■ - T n ;L)^. Basically, we 
consider a ID quantum conductor of length L, and com- 
pute the change of the transmission eigenvalues upon at- 
tachment of a thin slice of length 5L. Using the pertur- 
bation approach, we obtai n 11 ' 13 : 



1 2n 

- (ST,) = -V + N{N _ 1) 
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x \ l-T i + 2j2 Ti+ £-^ \ (5) 
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where 5s = SL/l, I is the mean free path defined by the 
first moment of the transmission eigenvalues of the thin 
slice: SL/l = l- (£? =1 T^SL)) /n. The third and higher 
moments vanish at the first order of Ss. DMPK equation 
is just the Fokker-Planker equation for the evolution of 
the distribution function P: 



8P X _ 2 " d I , . d fP x 
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where = TV(TV — l)/n, and we have re- 

expressed the distribution function in a new set of 
variables A* = (1 - T^/Ti, P A (Ai, A 2 • • • A n , s) = 
P(T U T 2 , ■ ■ ■ T n , L) JXi(l + A,)" 2 , and 
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which actually coincides with Eq. (j4]) except for an unim- 
portant denominator. We note that a similar DMPK 
equation was derived by Takane for the metallic carbon 
nanotubesii. 

Equations ([THE]) are the central result of this paper. 
The equation is reduced to the usual DMPK equation 
of the ordinary ID conductor of the symplectic ensemble 
(j3 = 4) for the even TV = 2n*±. On the other hand, for 
the odd TV = 2n + 1, the equation is modified, and one 
expects a different distribution of the transmission eigen- 
values. We will spell out its implications in the following. 

Equation ([7]) turns out to be completely integrable for 
both eveni£ and odd TV. To see this, we adopt a new set 
of variables {xi} that are related to {Xi} by A, = sinh 2 Xi, 
and P x ({xi}, s) = P\ 11"= 1 smn % x i- We further make the 
substitution P x = (x)^f(x, s), with 

71 

£(x) — YJ_ sinh 2 (2^ + Xj) sinh 2 (x j — Xj) 

Kj = l 

even TV 
sinh^ Xi odd TV 

and the equation is transformed to: 
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We can then identify the operator inside the square 
bracket of the rhs. of Eq. iflOj) being the radial part of 
the Laplace-Beltrami operator for the irreducible sym- 
metric space SO*(2JV)/J7(JV)lSi£. In particular, for the 
odd TV, Equation Q corresponds to a root system BC n 
of a e {±ei,±2ei 7 ±ei ± ej}, and has the appropri- 
ate multiplicity m a = {4,1,4} for a successful map- 
ping to the Laplace-Beltrami operator (see Table Bl of 
Ref. [lU l . This allows us to express the distribution func- 
tion P x (x,s) as a superposition of the zonal spherical 
functions &k(x)^: 



P x ({x},s) = C(s)£ 2 (x) J <i>k(x)e- k 



-fe 2 s/2 7N 



d n k 

Rfcf 



(11) 



For the odd TV, the Gindikin-Karpelevich formula 
(Eq. (C12) of Ref. El) yields: 
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where g(x) = ix(l + ix). The zonal spherical func- 
tion &k(x) can be constructed by a recurrent procedure 
(Eqs. (8.7-8.10) of Ref. M)- 

Using the asymptotic expansion of &k(x) and following 
the same line of derivations as shown in Ref. [HI, we can 



determine the asymptotic forms the distribution func- 
tion. In the localization regime (s ^> N): 

n 

P x ({x l },s)(xY[exp[-(j N /2s)(x l -x l ) 2 ] , (13) 
i=i 

where Xi = (s/7jv)[3+4(i— 1)]. It follows that the average 
conductance a ~ 1 + 2exp(— L/2£) with the localization 
length £ = Nl/2 for the odd TV, compared with £ = 
2(7V — 1)Z for the ordinary ID conductor (even AT)11. The 
conductance will have a log-normal distribution in this 
regime. 

In the diffusive regime (KsCiV): 



Px({xi}, s) cx TT (sinh 2 Xi — sinh 2 Xj) 2 (x 2 — x 2 ) 2 

x Yl e~ x ^ N/2s (x t sinh2a; i ) 1/2 (a* sinh a*) 2 . (14) 



Compared with the ordinary ID conducto r 15 ' 17 , 
the distribution function acquires an extra factor 
sinha^) 2 . Note that the correlations between the 
different transmission eigenvalues do not change. 

We can determine the average and variance of the con- 
ductance in the regime s <C TV using the method of mo- 
ments of Mello and Ston o 11 ' 18 , which computes the mo- 
ments of M q = J27=i as expansion in inverse powers 
of TV. From the DMPK equation ([?])■ we can establish a 
chain of the coupled evolution equations for moments of 
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K 

- 2pq (AfJMi) + pq 2 (MV- l M q+l ) - P q[qT 1] <M|) 
+ p(p - l)g 2 (MP- 2 (M 2q - M 2q+1 )) , (15) 

where =p sign stands for the odd (+) and even (— ) 
TV, respectively. Since MP = 0(N P ) in the particular 
regime, we can close the above equation order by order 
in the large TV limit. Noting that the average conduc- 
tance cr/uo = 2 (Mi) + Z 2 and the variance var(cr)/(7 2 = 



{(Ml) - (Afx) 2 ), cr = e 2 /h, we obtain: 
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+ 0(N~ 1 ), (16) 
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where Sa = a — Ncxq/{1 + s) is the weak localization cor- 
rection to the conductance, and we have introduced an 
index Z 2 that takes the value of (1) for the even (odd) 
N. Compared with the ordinary ID conductor, the edge 
transport of QSHI will have a different weak localization 
correction but the same universal conductance fluctua- 
tion. 




Figure 2: (color online) The weak localization correction 
to the conductance (Sa) and the variance of conductance 
(var(cr)). Left: The results obtained from the numerical solu- 
tions of DMPK equation for different values of iV (the num- 
bers by the curves). The N — > oo asymptotic behaviors 
Eqs. l|16H17p are indicated by the red curves. Right: The 
comparison between the results from the DMPK equation and 
that from the simulation of the 2D model Eq. I|18p . The re- 
sults for N = 4 and N = 5 are shown. In the 2D simulation, 
the parameters for the Kane-Mele model are: t = 1, A so = 0.2, 
X R = 0.2, A„ = (see Ref. [3 for the definitions of these pa- 
rameters), and the other parameters are: W = 0.6, ti = 0.2. 
The stripe has a width 48 lattice sites (LS) along the zig-zag 
direction. We have determined the mean free paths to be 
I = 2240 LS (N = 4) and I = 1680 LS (N = 5), respectively. 
200 disorder configurations are averaged. 



We have numerically solved the DMPK equation for 
different values of N, using a Monte-Carlo approach that 
simulates the diffusion of the transmission eigenvalues. 
The weak localization correction and the fluctuation of 
the conductance are calculated, shown in the left panel of 
Fig. [2 Equations lfT6|) - (fT7|l well predict the behaviors for 
both Sa and var(a) in the regime s <§; N. The difference 
between the odd N and the even N is evident. It is 
interesting to observe that although the variances of the 
conductance have the same asymptotic behaviors in the 
regime s <C N for the odd and even N (see Eq. (flT]) . 
they are very different in the crossover regime (s ~ N): 
var(cr) rapidly decreases to zero for the odd N, while it 
continues to increase and peaks at a value of 0.3(e 2 /ft.) 2 
for the even TV. 

We further test our results against a real model of the 
2D QSHI. We consider a system of a stack of TV lay- 
ers of honeycomb lattice, each of which is described by 
the tight-binding Hamiltonian introduced by Kane and 
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Figure 3: The distribution of conductance for N = 4 (lower) 
and 5 (upper). The drop lines shows the results from the 
DMPK equation. The solid lines indicates the distributions 
predicted by the asymptotic formulae Eq. l]14p [(a)-(c)] and 
Ip]) [(d)]. Atr = a(s) - tr(oo). 
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where cu (c' a ) denotes the annihilation (creation) opera- 
tor for lattice site i at the Z-th layer, iJ ; KM denotes the 
Hamiltonian for the Zth-layer and has the same form as 
Eq. (1) of Ref. 0, e u is the random site energy that uni- 
formly distributes in [— W/2, W/2], and the last term 
introduces the hopping between the neighboring layers 
with a hopping constant t\. With the appropriate pa- 
rameters, the system becomes an insulator with iV-pairs 
of edge states. The presence of the Rashba spin-orbit cou- 
pling in H KM and the disordered site energies will intro- 
duce backscattering between different edge channels. In 
the simulation, a small value of W is chosen, so that the 



bulk is still insulating. This is different from the previous 
numerical investigations which concern more on the an- 
nihilation of the edge states by the strong disorder due to 
the breakdown of the bulk gap 19 i 20 . We have adopted an 
iterative approach based on the non-equilibrium Green's 
function to calculate the conductance of a stripe of vary- 
ing length 2 ^. The results are presented in the right panel 
of Fig. [2] It is evident that both the weak localization 
correction to the conductance and the variance fit well 
with those predicted from the DMPK equation. It justi- 
fies our presumption that a 2D QSHI is effectively a ID 
quantum conductor, and can be described by the DMPK 
equation {2HHJ, 

Figure [3] shows the distributions of conductance for 
different values of s. A crossover from the Gaussian dis- 
tribution at the ballistic limit (s = 1) to the log-normal 
distribution in the localization regime (s = 10, 20) can be 
observed. The difference between the odd and the even 
N is the most notable in the crossover regime (s = 3, 5), 
where the distribution for N = 5 shows a smooth high 
conductance tail, while that for N = 4 has the a sharp 
threshold. 

Finally, we discuss the possible experimental verifica- 
tion of our theoretical predictions. A 2D QSHI with the 
multiple pairs of edge states can be realized by confining 
the 3D QSHI in one direction 22 . An alternative and more 
flexible way is to put an ordinary mesoscopic ID quantum 
wire in the proximity of the edge of a 2D QSHI sample, 
and couple them through a controllable gate. The en- 
semble statistics can then be measured for the combined 
system. By varying the coupling strength between the 
ID wire and 2D QSHI, one expects a crossover of the 
ensemble statistics from the even N to the odd N. 
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